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The stability and bonding of the ternary complex K2PtCl6 structure hydrides is discussed using 
first principles density functional calculations. The cohesion is dominated by ionic contributions, 
but ligand field effects are important, and are responsible for the 18-electron rule. Similarities to 
oxides are discussed in terms of the electronic structure. However, phonon calculations for Sr2RuH6 
also show differences, particularly in the polarizability of the RuHe octahedra. Nevertheless, the yet 
to be made compounds Pb2RuH6 and Be2FeH6 are possible ferroelectrics. The electronic structure 
and magnetic properties of the decomposition product, FeBe2 are reported. Implications of the 
results for H storage are discussed. 
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I. INTRODUCTION 

The complex hydrides, DMRq, D=Mg,Ca,Sr,Eu 
and M=Fe,Ru,Os form in the cubic (Fm3TO,#225), 
K2PtCl6 structure. ^~^ This structure has D on site 8c 
(1/4,1/4,1/4), M on 4a (0,0,0), and H on site 24e 
(a;H,0,0). These compounds are of fimdamental interest 
because of the unusual structural motif and the interest 
in understanding resulting electronic structure, and the 
bonding associated with it. Furthermore, they could be 
of practical interest as potential hydrogen storage mate- 
rials. Mg2FeH6 has very high volume and mass storage 
efficiency (150 g/1 and 5.4 wt. %), but is too stable for 
reversible room temperature applications.^'^'* In this re- 
gard, understanding of the electronic structure and cohe- 
sion may be helpful in finding modifications that improve 
the thermodynamics, to produce a material for hydrogen 
storage in mobile applications. 

The crystal structure may be regarded as a cubic 
double perovskite A2BB'}le, -with A = D, B=M and 
_B'=vacancy. Therefore, from a structural point of view, 
the compounds consist of A/He octahedra, well separated 
by presumably inert D ions, whose role is to fill space 
and donate charge to the A/Hg units. The cubic Fmim 
structure is maintained for all these compounds in spite 
of large variations of the A and B site cation radii, in 
contrast to the structural distortions often found in ox- 
ide perovskites and double perovskites. 

These compounds are insulators, and, like many of the 
complex hydrides, '^ follow the 18 electron rule, which says 
that the number of non-bonded metal electrons plus the 
number of electrons in the metal-ligand bonds should be 



18. In the simplest view, this would correspond to full 
s, p and d shells associated with the [AfHg]"'" structural 
units. However, Miller and co-workers,^ have emphasized 
the importance of ligand field effects in these complex 
hydrides, and, in fact, the calculated electronic struc- 
tures for these materials show band gaps within the d 
manifolds, indicating a more complex situation. ^^ In par- 
ticular, insulating band structures, in qualitative accord 
with experiment, resulting from band gaps between crys- 
tal field split transition metal d manifolds were found in 
non-self-consistent warped mufhn-tin calculations, based 
on the Xa method with a = 1. 

Here, we re-examine the electronic structures, which 
we obtain using a full-potential, self-consistent linearized 
augmented planewave (LAPW) method, and use these 
results, along with calculations of the formation en- 
thalpies, to discuss the bonding of these materials and 
possible directions for modifying them to alter their sta- 
bility. In addition we discuss the hypothetical materials, 
Pb2RuH6, which was studied as a potential ferroelectric, 
and Be2FeH6, which is both in relation to ferroelectricity 
and to better understand their stability. 



II. APPROACH 

As mentioned, the calculations were done within 
the local density approximation (LDA) using the gen- 
eral potential linearized augmented planewave (LAPW) 
method. ^^ Local orbital extensions'^ were used to re- 
lax the linearization errors for the transition metal d 
states, and to treat the upper semicore levels of the 



alkaline earth and transition metal atoms. An LAPW 
sphere radius th = 1-1 ao was used for H in all the com- 
pounds. Metal radii of 1.8 ao were used for Mg2FeHe 
and Ca2FeH6. For Sr2FeH6, metal radii, rsr = 2.0 oo 
and rpc = 1-8 flg were used for Sr and Fe, respectively. 
For Ca2RuH6 and Sr2RuH6, metal radii of 1.95 aq were 
used. For Mg2RuH6 we used rug = 1-8 ao and tru = 1.9 
ao. Well converged basis sets, defined by raGmax—Q-O, 
where Gmax is the planewave cut-off were employed. 
The Brillouin zone sampling during the iteration to self- 
consistency was done using the spccial-k points method 
with a 4 X 4 X 4 mesh, which corresponds to 10 k-points in 
the irreducible wedge. Densities of states were generated 
using a tetrahedron mesh of 145 k points in the wedge. 
Convergence was tested, both for the zone sampling and 
basis set size, by repeating some calculations with higher 
raGmax^T and more k points (8x8x8). Based on these 
tests, the present convergence with respect to these pa- 
rameters is better than 2 mRy/cell for total energies, and 
better than 1 niRy for band energies. 



III. ELECTRONIC STRUCTURE AND THE 
EIGHTEEN ELECTRON RULE 



Crystal field effects, which are important here, may be 
expected to be sensitive to the H positions. Here we use 
LDA relaxation to determine the H positions (given in 
Table I). As may be seen, they are in good agreement 
with experiment for those compounds for which neutron 
refinements are available. Table I also gives the full sym- 
metry Ag Raman phonon frequency associated with the 
H internal structural parameter. Raman and inelastic 
neutron scattering measurements^'* for Mg2FeHg yield an 
experimental frequency of 1873 cm^^, in good agreement 
with the present LDA value of 1923 cm'^. The LDA 
frequencies follow the reported trend for H bond stretch- 
ing modes in infrared data,^ decreasing with increasing 
D ionic radius, and increasing with increasing M atomic 
number. It is interesting to note the increase in frequency 
from Sr2RuHe to Sr20sHe. Generally, the ionic proper- 
ties of Ad and 5d elements are very similar due to relativis- 
tic contraction. The main difference, which is also due 
to rclativistic contraction, is that the s states are lower 
relative to the d states in the 5d series. This leads to a 
higher position of the d states relative to ligand states 
in compounds where the transition elements are cations. 
The result is reduced covalency. Since covalency softens 
ionic interactions, the result is a stiffer lattice in 5d ionic 
compounds relative to the corresponding 4d compounds. 
A good example is the comparison of KNbOa , which is a 
good ferroelectric with KTaOs, which has practically the 
same lattice parameter, but is not ferroelectric. *° 

The calculated electronic band structures are shown in 
Figs. 1, 2 and 3. The corresponding electronic densities 
of states (DOS) are in Figs. 4, 5 and 6. The electronic 
structures are qualititively similar to those of Orgaz and 



Gupta (Ref. 11), in that all the materials are insulating 
with band gaps within the transition metal d bands. 

The band structure for all the compounds studied con- 
sists of 6 H Is derived bands, holding 12 electrons, fol- 
lowed by crystal field split transition metal d bands. The 
Is band width is largest in the Mg compounds, and in 
the case of Mg2RuHe this is sufficient to yield an overlap 
between the Is manifold and the lower d manifold. In the 
octahedral crystal field of the H, the metal d bands sep- 
arate into a 3-fold degenerate (6 electrons) t2g manifold 
and a 2-fold degenerate (4 electrons) eg manifold. The 18 
valence electrons then populate the H derived Is bands 
and the metal t2g bands; the insulating gap is between 
the occupied t2g and the unoccupied Cg manifolds. Since 
these are gaps within a crystal field split d-band, LDA 
band gap errors are expected to be small in the absence 
of strong correlation effects.*^ However, we note that the 
gap is between narrow d bands, and so a larger gap may 
be observed both because of optical dipole selection rules 
and, in the Fe compounds, correlation eflfects. 

The His character of the lowest 6 bands may be seen 
from the projections of the DOS onto the LAPW spheres. 
The small 1.1 ao H spheres, used here, imply that most 
of the charge of H ions will be outside the sphere. A 
free H~ ion, stabilized by a Watson sphere (to approx- 
imately represent the Madelung field) of radius 3.01 A, 
would have only 0.45 e (of 2) inside a sphere of radius 
1.1 ao. Integrating the H s projection of the DOS over 
the lowest 6 bands, we obtain from 0.50 e (Ca2RuHg and 
Sr2RuH6) to 0.52 e (Sr2FeH6) inside each H sphere, not 
far from this simple ionic view especially if one allows for 
a somewhat different breathing. This is also similar to 
what was found for NaAlH4.-'^^ Thus the basic electronic 
structure is ionic consisting of H anions and D and M 
cations. These compounds should therefore be viewed as 
ionic for understanding the crystal cohesion. 

However, from the stand-point of understanding the 
electronic structure and H storage properties, covalency 
is important. The effects of M - H hybridization are 
clearly seen in the electronic structures. While the bot- 
tom 6 bands are essentially H Is bands, they contain 
the two formally bonding H s - M Cg a combinations 
(as well as four non-bonding combinations). The lowest 
band is the symmetric combination of s orbitals, which 
is a non-bonding combination with M Cg states, but is 
formally bonding with the nominally unoccuppied M s 
states. This is followed by five more H s bands, including 
the non-bonding and the formally bonding H s - M e^ 
(7 combinations. In terms of degeneracies, this division 
into a "s"-like one fold symmetric band, a five- fold man- 
ifold and the three-fold t2g manifold is formally like the 
s, p, d electron counting of the 18-clcctron rule. How- 
ever, this counting does not correspond to atomic level 
filling. Instead, the 2-electron "s" and 10-electron "rf" 
manifolds are actually from combinations of hydrogen s 
states, and are therefore very different in atomic char- 
acter from the 6-electron, "p" group, which comes from 
very weakly hybridized M t2g bands. 



The lowest conduction bands derive mainly from the 
corresponding antibonding combinations and sd derived 
states associated with the D cation. The covalency in 
the [MHg]'*^ units can be seen in the M d contributions 
to the H s bands, for example. The i2g bands show much 
less hybridization, as expected in an octahedral ligand 
environment. It should be emphasized that the bands 
show relatively little dispersion, with the exception of the 
D=Mg compounds, and that there are generally clean 
gaps between the different manifolds (H s, t2g and Cg), 
which implies that weakly interacting [MHg]^" units may 
be regarded as the basic building blocks for understand- 
ing the band structure. The sizable crystal field splitting 
of the M d bands underlies the 18-electron rule in these 
compounds. In particular, without it the t2g and Cg man- 
ifolds would overlap, and then there would be no barrier 
to adding more than 18 electrons; the 18 electron rule 
here is the result of the crystal field splitting of the metal 
d levels. The octahedral geometry, with its large ligand 
field is energetically favorable for 18 or fewer electrons. ^^ 
The substantial crystal field splittings (as compared to 
the on-site Hund's coupling, which can be characterized 
by a Stoner / ~ 0.7 — 0.9 eV for Fe)^*'^^ are responsible 
for the low spin Fe observed in these compounds. 

To summarize the results so far, the electronic struc- 
ture is built up in the following way in decreasing order of 
the size of the interactions involved. (1) Coulomb inter- 
actions, particularly the Madelung field, stabilize an ionic 
configuration, nominally D^^ M^^Hg . This is the main 
ingredient in the cohesion. (2) Hybridization between 
the H s orbitals and the M Cg orbitals lead to a bonding 
anti-bonding splitting between these and contribute to 
a substantial crystal field splitting between weakly hy- 
bridized occupied M <2g states and unnoccupied M eg 
states. This splitting and the position of the D derived 
states well above the t2g energy underlies the 18 electron 
rule. (3) Hopping between the [Mng]**" units (presum- 
ably mostly assisted hopping via unnoccipied D s and d 
states) leads to band formation. This is reminiscent of 
some of the oxide double perovskites, A2MM'0e, with 
an inert M', such as Sr2RuY06, although in that case 
the hybridization inside the [RuOg]''^ units is very much 
stronger than in the present hydrides. ^° 



IV. PHONONS, FERROELECTRICITY AND 
HYPOTHETICAL PbaRuHe 

The resulting picture of ionic crystal with substantial 
covalency between anions and an octahedrally coordi- 
nated transition element cation suggests similarities with 
double perovskite oxides. Furthermore, the fact that the 
latttice contains a large anion (H~) stabilized by the 
Madelung field and hybridized with nominally unoccu- 
pied transition metal states further suggests connections 
with perovskite oxides, particularly ferroclectrics. In 
fact, many of the technologically important ferroelectrics 
are based on perovskites ABO3 with mixtures of metal 



atoms on the B sites (these can be disordered or ordered 
as in e.g. double perovskite). Examples include PZT 
[Pb(Zr,Ti)03], PMN-PT [Pb(Mg,Nb,Ti)03] and PZN- 
PT [Pb(Zn,Nb,Ti)03]. In these materials both Pb-0 and 
B-0 hybridization is important in the ferroelcctricity.'^^ 

In order to further elucidate the relationship to oxides, 
we calculated those zone center phonon frequencies of 
Sr2RuIIg compatible with a rhombohedral i?32 symme- 
try, and compare with similar calculations for hypothet- 
ical Pb2RuH6. This was done at the experimental lat- 
tice parameter of Sr2RuHg, using the relaxed H position. 
This non-centrosymmetric group would include the fer- 
roelectric mode, if the material were ferroelectric. The 
calculations were done by fitting the dynamical matrix 
to a series of frozen phonon calculations with small dis- 
placements of the various atoms. This yields 6 three-fold 
degenerate modes (plus the three w = acoustic modes) . 
The frequencies and displacement patterns of the phonon 
modes are given in Table IV. The 2 % difference between 
the Ag frequencies for Sr2RuH6 between Tables I and IV 
refiect the different approaches and should be considered 
indicative of the errors in the fits used in constructing the 
dynamical matrix. The fitting errors can also be seen in 
the deviation of the mode character as given in Table IV 
from the Ag character required by symmetry. For the 
Ag frequency, the value in Table I should be considered 
more reliable because that value was obtained enforcing 
the exact mode symmetry, but it should be kept in mind 
that the LDA error is likely larger than the difference 
between the values in Tables I and IV. 

The highest frequency branches correspond to H mo- 
tions, as expected. Of these, the Ag Raman mode, which 
corresponds to symmetric breathing of the RuHg octahe- 
dra, is the stiffest mode, and the next lower mode also 
involves modulation of the Ru-H bond lengths. The two 
intermediate modes (738 cm^^ and 777 cm^^) involve 
distortion of the octahedra, which would also yield lower 
frequency modes in oxides. The two lowest frequency 
modes are motions of the Sr atoms within their cages. 
The lowest mode is the antisymmetric motion of the Sr, 
which is not ferroelectric. The frequencies of these Sr mo- 
tions are compatible with the frequencies of the shearing 
modes that modulate Sr-H distances when the mass dif- 
ference is accounted for. 

This pattern of phonon modes is quite different from 
what would occur in an oxide near ferroelectricity. In 
that case, there would be a low frequency cooperative 
mode. This would consist of a distortion where the 
cations move relative to the O atoms comprising the oc- 
tahedra, reflecting the high polarizability of the octahe- 
dra softened by covalent interactions. '^^''^^ Here, the mode 
corresponding to motion of the Ru with respect to the H 
is at high frequency (1433 cm^^) and the lower frequency 
Sr derived modes have only a small component of Ru mo- 
tion relative to the H octahedra. 

We repeated the calculations for the hypothetical com- 
pound Pb2RuH6. In perovskite and double perovskite 



oxides, Pb typically can be substituted for Sr. The Pb 
compounds typically have unit cell volumes very close to 
the Sr analagues, but are more likely to be ferroelectric 
because of Pb-0 covalency (e.g. PbTiOa vs. SrTiOs). 
Since the experimental lattice parameter of Pb2RuHe is 
unavailable, we used the value for Sr2RuHe. This choice 
also makes comparison of the two systems more direct. 
LDA relaxation of the lattice parameter yielded a value 
1.8% smaller than this, but considering the usual under- 
estimate of lattice parameters in the LDA, we do not 
consider the LDA value to be more reliable than the use 
of the SraRuHg valuers The PbaRuHg modes (Table V) 
are qualitatively like those of Sr2RuH6, except that the 
low frequency Pb modes are shifted down in frequency. 
This shift is, however, larger than can be accounted for 
by the mass difference, and the modes are reversed in or- 
der. The lowest mode is now the symmetric mode and it 
is at zero frequency to within the precision of the present 
calculations. We also calculated the energy as a function 
of the rotation of the RuHg octahedron. Such rotational 
degrees of freedom compete with ferroelectricity in fer- 
roelectrics such as Pb(Zr,Ti)03.^'* Here, these modes are 
stable. We obtain frequencies of 576 cm~^ and 490 cm~^, 
for Sr2RuH6 and Pb2RuH6, respectively. This suggests 
Pb2RuH6 as a candidate ferroelectric hydride. We found 
a similar result for hypothetical Be2FeH6. In this case, 
we obtained a slightly unstable mode of Be character, 
with a ferroelectric displacement pattern, and stable ro- 
tational and antiferroelectric modes. If these compounds 
are made, the possibility of ferroelectricity should be in- 
vestigated, e.g. by low temperature structural studies 
and temperature dependent dielectric measurements. 



V. ENERGETICS AND ZERO POINT EFFECTS 



In order to better understand the stability of these 
compounds, we performed calculations of the formation 
enthalpies by comparison of the total energies with those 
of decomposition products. Specifically, we did calcula- 
tions for the elements, Mg,Sr,Ca,Fe,Ru and Os in their 
bulk metallic form (in the LDA at the experimental lat- 
tice parameters, parallel to the calculations done for 
the hydrides, including ferromagnctism for Fe), the H2 
molecule (relaxed, in the LDA) and MgH2, Call2 and 
SrH2 (using experimental lattice constants, but relaxed 
atomic positions) . In addition calculations were done for 
elemental Be and the intermetallic phase Be2Fe, which is 
the expected decomposition product of the hypothetical 
phase Be2FeH6 (see below). 

For the H2 molecule we used a cubic supercell of lattice 
parameter 4.5 A. This yielded an LDA energy of -2.288 
Ry, and bond length of 0.765 A, and bond stretching 
vibrational frequency of 4217 cm^^. These results are 
in good agreement with previous LDA calculations. For 
example, Patton and co-workers report a vibrational fre- 
quency of 4188 cm~i and bond length of 0.765 A.^^ The 



H2 zero point energy obtained from the LDA frequency 
is 25.2 kJ/mol. This is a substantial number, which un- 
derscores the fact^^^^^ that zero point effects need to be 
considered in the thermodynamics of hydride formation. 
Here we neglect metal modes, and consider only the H 
contribution to the zero point energies, which we write 
as ihuj per H2 unit, where Co is an average H vibrational 
frequency. The effective uj for H2 is 703 cm~^, so for hy- 
drides with w > 703 cm~^, zero point motion will reduce 
the formation energy, and the corresponding deuterides 
and tritides will form more easily than the hydrides, while 
for materials with oj < 703 cm~^, the converse will be 
true.^^ At least in principle, this difference can be used 
to obtain the average H frequency from experimental for- 
mation energies of hydrides and deuterides, but to our 
knowledge this has not been done for these materials. 

In order to estimate lu for the compounds considered 
here, we performed LDA calculations for selected distor- 
tions and assumed that the H behaves in an Einstein-like 
way. In particular, for MgH2, we displaced a single H in 
the unit cell (which contains 4 equivalent H atoms) along 
the three principal directions in its cage and averaged 
the resulting frequencies to obtain an average frequency 
(the principal directions relative to the lattice are [1,1,0], 
[1,1,0] and [0,0,1]). For CaH2 and SrH2, we used the av- 
erage of the four full symmetry Raman modes that are H 
derived (there are two other such modes associated with 
metal motion). For the K2PtCl6 structure hydrides, we 
used the average of the highest two F-point modes consis- 
tent with rhombohedral symmetry (the highest of these 
is the full symmetry Raman breathing mode) to obtain 
an effective B-site metal - H bond stretching force con- 
stant (which contributes 1/3 of the modes) and obtained 
the effective force constant for the other two thirds of the 
modes by averaging the two other H derived modes from 
the rhombohedral symmetry and the octahedral rotation 
mode. As a test, we also calculated an average Einstein 
frequency for Mg2FeH6 by displacing a single H (of the six 
in the unit cell) perpendicular to the Fe-H "bond", and 
along it, similar to the procedure that was followed for 
MgH2. This yielded a shear frequency of 727 cm^^ and 
a stretch frequency of 1828 cm^^, for an average tt'=1094 
cm""'^, in fortuitously good agreement with the estimate 
of 1089 cm^^, made as above (Table HI). In all cases, 
the averages are arithmatical averages of frequencies as 
is appropriate for the zero point energy. 

a-MgH2 has a tetragonal structure (spacegroup 
P4:2/mmm) with one H coordinate, x. We obtain 
a;=0.3046, in agreecment with the recent neutron mea- 
surement of Bortz and co-workers, ^^ who obtained 
a;=0.3040, and an LDA calculation by Yu and Lam,^^ 
who also obtained a;=0.304. Moving a single hydro- 
gen in the unit cell of MgH2, we obtained frequencies 
of 1277 cm~^ along [1,1,0] (bond stretching) 592 cm~^ 
along [1,1,0] (bond bending) and 993 cm"-'^ along [0,0,1] 
(mixed). The bond stretching "Einstein" frequency is 
near that of the full symmetry Raman mode, which we 
obtain at 1301 cm^^ and which is also of bond stretch- 



ing character. This supports the simple Einstein-hke ap- 
proach used. 

Good agreement with neutron diffraction results^^'^^ 
is also obtained for the internal structural parameters of 
CaH2 and SrH2, which occur in an orthorhombic Pnma 
structure, as given in Table II. As mentioned, the full 
symmetry Raman frequencies obtained from this relax- 
ation, were used to construct the H frequency for the zero 
point contribution to the enthalpy of these compounds. 

LDA formation energies are given in Table III. The for- 
mation energies of MgH2, CaH2 and Srll2 are in excellent 
agreement with experiment and also in good agreement 
with previous LDA calculations.^^ The formation energy 
of Mg2FeIl6 is the best studied of the K2PtCl6 hydrides, 
and is reported as -98 kJ/mol (Ref. 2), -86 kJ/mol (Ref. 
7), and -77.4 kJ/mol (Ref. 8). The calculated LDA en- 
ergy of -133 kJ/mol is significantly larger and this dif- 
ference would seem to be at the high end of the nor- 
mal range of LDA errors, especially considering the good 
agreement with experiment for Mgll2, and FeAl, where 
we obtain agreement with experiment to within 5 kJ/mol 
of Fe.^° Besides LDA errors, the most likely source of er- 
ror is the crude method that we used to obtain an average 
phonon frequency. However, even if the average phonon 
frequency were 250 cm~^ higher than our estimate, which 
we think is unlikely, the calculated enthalpy would shift 
by only 9 kJ/mol of II2. 

One possibility is that some of the difference is experi- 
mental in origin, related to the possible existence of some 
stable hydride among the decomposition products, which 
would then stabilize the products relative to Mg2FeIl6 
and therefore would lower the formation energy as mea- 
sured by the decomposition. In any case, we do find 
certain trends. First of all, Mg2FeIl6 is by far the least 
stable of the K2PtCl6 hydrides studied. However, this 
is connected with the fact that MgH2 is much less sta- 
ble than SrH2 or CaH2. If one considers formation via 
2£)H2 + II2 + M —> 1)2 MHe, then this heat of forma- 
tion is largest for Mg2FeH6 as might be expected from 
an ionic picture. Secondly, the formation energy per H2 
is significantly larger for Mg2FeH6 than for MgH2. This 
implies that under thermodynamic conditions, without 
some very unusual entropy contribution, the decomposi- 
tion should proceed directly to H2 and the elements, or, 
at a bare minimum, if there is an intermediate hydride 
phase, it should not be MgH2. The formation energies 
of the Ca and Sr compounds on the other hand are very 
close to those of CaH2 and SrH2, so depending on the 
conditions, those decompositions may very well proceed 
via an intermediate M + DH2. 

Similar calculations were done for the other com- 
pounds in order to estimate the H zero point energy, but 
these were at a lower level of convergence in the fitting 
of the dynamical matrix. 



VI. STABILITY, BONDING AND IMPLICATIONS 

We now speculate about possible implications of our 
results for hydrogen storage. First of all, we note that the 
cohesion is ionic, and that it is the Madelung field that 
stabilizes the [MHg]^" units. Changes in the Coulomb 
potential then ought to strongly affect the bond lengths 
in these units as well as their stability. This is already 
apparent in the values of a;H of Table I, which show sub- 
stantial changes in the M-H bond lengths as the lattice 
parameter is changed by D site substitution (note that 
the octahedra are not connected so that they need not 
breath with the lattice) . This implies significant tunabil- 
ity in the properties with substitutions. Secondly, the 
ionic stabilization of the lattice implies that mixed M 
cation substitutions should be possible and that the octa- 
hedral coordination of the metal atoms will be preserved 
in them. For example, if the partial or full substution Os 
^ ^ Re -I- ^ Ir could be made, its structure is expected to 
feature ReHg and IrHg octahedra, rather than different 
Re-H and Ir-H coordinations on this lattice. However, it 
is unclear if any of these substitutions can be made, and 
even if they can it is unclear whether they will be benefi- 
cial. Finally, we note that the fully hydrided compound 
can therefore be destabilized by driving the transition 
metal d states up in energy via the Madelung potential, 
if all other things were equal. One way would be to sub- 
stitute some fraction of the Mg with a monovalent cation 
if one can be made to enter the lattice. 

However, the stability is relative to the decomposed 
products, and it is clear from the calculated energetics 
that these play a major role. For example, Mg2RuHg 
and Mg2FeH6 are the least stable compounds relative to 
decomposition into elemental Mg and Ru or Fe, but they 
are the most stable with respect to a hypothetical inter- 
mediate MgH2 + Fe/Ru. As mentioned, this suggests 
that under normal conditions Mg2FeH6 and Mg2RuH6 
decompose into H2 and elemental metals without any 
MgH2 intermediate, consistent with the observation of 
Bogdanovic and co-workers.* 



VII. REDUCING THE FORMATION ENERGY: 
HYPOTHETICAL BeaFeHg 

Considering the trend in the energetics with respect to 
the alkaline earth element, one possibility for obtaining a 
lower formation energy would seem to be replacement of 
Mg by Be. This would seem especially likely considering 
the properties of Be metal, which include strong bond- 
ing that would compete with the formation of hydride 
phases. In order to check this trend we performed calcu- 
lations for hypothetical Be2FeH6 to obtain its formation 
energy. Since the compound is hypothetical, we obtained 
the lattice parameter by relaxation in the LDA. The cal- 
culated structure has a lattice parameter of 5.65 A, a 
H internal coordinate of x=0.2648, and a corresponding 



full symmetry Raman phonon frequency of 2233 cm"-'^. 
No doubt the LDA underestimates the lattice parameter 
slightly, as is typical. In any case, with the LDA struc- 
ture, we obtain a static formation enthalpy of -37 kJ/mol 
H2 with respect to elemental products. This confirms 
the conjecture that Be would lead to lower binding ener- 
gies. However, while this energy suggests that Be2FeH6 
would be an interesting hydride phase, it neglects the fact 
that unlike Mg, Be forms compounds with Fe. In par- 
ticular, FeBe2 is a known intermetallic compound and 
would compete with the hydride phase. We calculated 
the enthalpy of formation of FeBe2 (details are in the 
next section) and find -87 kJ/mol. The relevant energy 
for the stability of the hydride Be2FeH6 is therefore not 
the formation enthalpy from the elements, but from the 
intermetallic, BeFc2, which, with H2, would be the prod- 
uct of the decomposition. Relative to decomposition into 
the elements, the existence of the intermetallic then re- 
sults in a shift of the formation enthalpy of Be2FeHg by 
29 kJ/mol on a per H2 basis, to yield -8 kJ/mol. The 
average phonon frequency, determined as for the other 
K2PtCle hydrides, discussed above, is 1172 cm~^, yield- 
ing a zero point correction of -|-17kJ/mol H2, placing the 
calculated enthalpy including zero point at +9 kJ/mol. 
Thus, it is likely that Be2FeH6 is marginally unstable 
with respect to decomposition into FeBe2, and therefore 
will only be formed under high pressure or by chemical 
routes. 



VIII. ELECTRONIC STRUCTURE AND 
MAGNETISM IN FeBez 

FeBe2 is an interesting hard magnetic material. In 
particular, it has a relatively low density, high anisotropy 
and a very high Curie temperature, Tq of 880 K.^^~'^^ Ex- 
perimentally, FeBe2 crystallizes in the hexagonal MgZn2 
structure (spacegroup PG^/mmc, No. 194) and has mag- 
netization, M « 1.95^B/Fe. Since we are not aware of 
previous first principles studies of this material, we briefly 
summarize our results for the electronic structure of this 
compound. 

The reported lattice parameters are a=4.215 A, and 
c=6.853 A.^^ The unit cell contains four formula units. 
The Be atoms are on sites 2a (0,0,0) and 6h (x,2x,l/4), 
while the Fe atoms are on 4/ (1/3, 2/3, z). Experimental 
values of the two internal parameters are not available 
from experiment, so they were found by structural re- 
laxation in the LDA. We find x= 0.8294 and z=0.061. 
The calculated LDA spin magnetization is 1.76 /is/Fe 
and the magnetic energy is 0.237 eV/Fe. This is only 
~ 3kTc suggesting some itinerant character. 

The local spin density approximation (LSDA) band 
structure and density of states are given in Figures 7 and 
8. The band structures shows narrow crystal field split 
Fe 3d bands on top of a broad manifold of free electron 
like Be sp derived bands. These Be derived bands are 



weakly polarized, opposite to the Fe polarization, simi- 
lar to the case of YFe2, for example. ^^'^^ The majority 
spin Fe 3d bands are fully occupied, while the Fermi en- 
ergy falls in the crystal field gap between the i2g and eg 
manifolds in the minority spin. This yields two minority 
spin Cg holes per Fe, and explains the ^ 2 fis/Fe mag- 
netization. Relative to bcc Fe, there is a transfer of Fe 
s character to the Be derived bands, and a back trans- 
fer of charge to give effectively neutral Fe, with eight 3d 
electrons. This is consistent with the picture discussed 
by Jesser and Vincze based on experimental susceptibil- 
ity and Mossbauer measurements.^® This pseudogapped 
band structure yields a relatively low density of states 
at the Fermi energy in both spin channels, N^{Ep)=0.39 
eV^^ and Ni{Ep)=0.67 eV~^ on a per formula unit ba- 
sis. 

The calculated formation energy of FeBe2 is -87 kJ/mol 
on a per formula unit basis. As a test of our approach, 
we also calculated the formation energy of FeAl in the 
same way. The result was -77 kJ/mol(FeAl), which is 
in good agreement with the experimental value of -72.6 
kJ/mol.'^"'"'" This suggests that the error in the formation 
energy of FeBe2 is likely in the range of 5 kJ/mol. 



IX. SUMMARY AND DISCUSSION 



The present density functional calculations show that 
the K2PtCl6 hydrides are ionic compounds, with some 
covalency. The 18 electron rule is a consequence of lig- 
and field effects on the transition metal site. This type 
of ionic character suggests the possibility of ferroelectric- 
ity in related hydrides. We find that the hypothetical 
compounds, Pb2RuH6 and Be2FeH6 are on the border- 
line of ferroelectricity, and should be investigated in this 
context, if they can be made. Further, the ionic char- 
acter stabilizes H anions, and Fe cations, which is why 
Fe participates in hydride formation, although Fe metal 
is not a hydride former. The ionic character implies a 
certain degree of tunability of the properties of these hy- 
drides, which may allow adjustment of the thermody- 
namics. However, since Fe and Mg do not form inter- 
metallic compounds, it is likely that the properties of 
Mg2FeH6 cannot be made better for vehicular H storage 
than those of MgH2, since MgH2 will be a competing 
phase. 

One way to reduce the stability of the hydride without 
facing this limitation would be to focus on the stability 
of decomposition products. One possibility would be to 
explore minor additions, X that are soluble in and stabi- 
lize an Fe-Mg-X intermetallic. These need not enter the 
hydride lattice, provided that they are available e.g. on 
hydride particle surfaces to promote the decomposition 
and provide sufhcient enthalpy via the formation of the 
intermetallic. This may be the most promising avenue for 
modifying Mg2FeH6 for hydrogen storage applications. 
As far as we are aware, the solubility of Mg in FeBe2 is 



not known. However, if Fe-Be-Mg intermetallics are sta- 
ble, the present results suggest that the addition of Be 
to Mg2FeH6 may lead to lower stability, which if not for 
the toxicity of Be, would be favorable for applications. 
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TABLE I. LDA and experimental H positions (a;H(LDA) 
and a;H(EXP), respectively), fully symmetric Ag Raman fre- 
quency, CO, and band gap. Eg. 



G. Tanaka, T. 
Alloys Compd. 



Mit- 
302, 





a;H(LDA) 


xh(EXP) 


a;(cm "*") 


Eg{eY) 


Mg2FeH6 


0.2412 


0.2420 [ 2 ] 


1923 


1.73 X-X 


Ca2FeH6 


0.2257 


0.2300 [ 4 ] 


1759 


1.27 X-X 


Sr2FeH6 


0.2178 




1694 


1.09 r-x 


Mg2RuH6 


0.2527 


0.2524 [ 4 ] 


1977 


2.93 X-X 


Ca2RuH6 


0.2359 




1837 


2.29 X-X 



SraRuHe 
SraOsHe 



0.2254 
0.2262 



0.223 [ 1 ] 



1778 
1893 



2.06 F-X 
2.26 r-x 



TABLE II. LDA and experimental atomic positions 
(a;(LDA), z(LDA) and a; (EXP), z(EXP) for Pnma CaHj and 
Srll2. j/=l/4 for all atoms in this structure. 

a;(LDA) 2(LDA) x(EXP) 2(EXP) 



CaHa Ca 


0.2380 


0.1100 


0.2378 


0.1071 


■271 


CaH2 HI 


0.3566 


0.4274 


0.3573 


0.4269 


[271 


CaHa H2 


0.9741 


0.6773 


0.9737 


0.6766 


[271 


SrHa Sr 


0.2382 


0.1109 


0.2438 


0.1108 


[281 


SrHa HI 


0.3558 


0.4278 


0.3570 


0.4281 


[281 


SrHa H2 


0.9732 


0.6787 


0.9693 


0.6825 


[281 



TABLE HI. LDA energies of formation on a per H2 basis, 
assuming full decomposition into separated elemental metals 
and H2- (minus means that the formation is exothermic). 
AHatatic denotes the LDA energy with no correction for zero 
point motion, uj is the average phonon frequency estimated 
from LDA calculations (see text) and A_H'(kJ mol~^) is the 
zero-point corrected formation energy for the hydride. The 
LDA vibrational frequency of H2 is used in this calculation. 



Compound 


AHatau.{^3 mol-1) 


Lu{cm """) 


AH{kJ mor^) 


Mg2FeH6 


-147 


1089 


-133 


CaaFeHe 


-221 


1033 


-209 


SraFeHe 


-211 


871 


-205 


Mg2RuH6 


-147 


968 


-137 


CaaRuHs 


-229 


1013 


-218 


Sr2RuH6 


-221 


1006 


-210 


Sr20sH6 


-216 


1023 


-205 


MgH2 


-90 


954 


-81 


CaHa 


-219 


767 


-217 


SrH2 


-207 


670 


-208 



TABLE V. Phonon frequencies and displacement patterns, 
as in Table IV but for PbaRuHe, a = 14.361ao; frequencies. 



cj are in cm 



Qpb 



A 



Pb 



7Ru 



eni 



eH2 



23 

52 
552 
632 
1265 
1626 



-0.023 
0.049 
-0.002 
0.002 
0.000 
-0.000 



-0.019 
-0.048 
-0.002 
-0.002 
0.000 
0.000 



0.089 
0.005 
-0.010 
-0.001 
0.013 
0.002 



0.082 
0.005 
-0.062 
-0.022 
-0.646 
-0.749 



0.038 
0.025 
0.518 
-0.473 
-0.028 
0.002 



0.077 
0.004 
-0.013 
0.041 
-0.745 
0.654 



0.045 
-0.017 
0.469 
0.520 
-0.008 
-0.040 



TABLE IV. Calculated frequences uj (cm^^) and dis- 
placement patterns for zone-center modes of Sr2RuH6, 
a — 14.361ao compatible with _R32 symmetry. The displace- 
ments are {asr,asr,ctSr) for the first Sr (in i?32 symme- 
try), {I3sr,l3sr,l3sr) for the second Sr, {'yRu,-fRu,'yRu) for 
Ru, {Shi,^hi,^hi) for the first type (in i?32 symmetry) of 
H (three atoms), and {5h2,^h2,^h2) for the second type of 
H (three atoms). The 5h are Ru-H bond stretch coordinates, 
while en are Ru-H bond shears. 



UJ 


QSr 


/3s. 


7Ru 


(5hi 


EHl 


(5h2 


eH2 


126 


0.073 


-0.078 


0.002 


0.001 


0.013 


0.002 


-0.009 


141 


-0.051 


-0.045 


0.075 


0.075 


0.037 


0.071 


0.041 


738 


0.003 


-0.000 


-0.005 


-0.024 


-0.673 


-0.079 


0.195 


777 


-0.002 


-0.003 


-0.011 


-0.087 


0.188 


-0.012 


0.670 


1433 


0.000 


0.000 


0.015 


-0.632 


-0.055 


-0.750 


-0.020 


1742 


0.000 


0.000 


0.002 


-0.757 


0.012 


0.643 


-0.044 
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FIG. 1. Band structure of Mg2FeH6 (top), Ca2FeH6 (mid- FIG. 2. Band structure of Mg2RuH6 (top), Ca2RuH6 (mid- 

dle) and Sr2FeH6 (bottom), with the relaxed LDA H posi- die) and Sr2RuH6 (bottom), with the relaxed LDA H posi- 
tions, tions. 




FIG. 3. Band structure of Sr20sH6 with the relaxed LDA 
H positions. 
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FIG. 4. (color online) Electronic density of states (heavy 
red line) and projections onto the LAPW spheres of Fe d char- 
acter (dashed green line) and H s character (dotted blue line) 
for Mg2FeH6 (top),Ca2FeIl6 (middle) and Sr2FeH6 (bottom), 
with the relaxed LDA H positions, on a per formula unit both 
spins basis. 
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FIG. 5. (color online) Electronic density of states (heavy 
red line) and projections onto the LAPW spheres of Ru d 
character (dashed green line) and H s character (dotted blue 
line) for Mg2RuH6 (top), Ca2RuH6 (middle), and Sr2RuH6 
(bottom), with the relaxed LDA H positions, on a per formula 
unit both spins basis. 





FIG. 7. (color online) Band structure of ferromagnetic 
FeBe2. Majority (minority) spin bands are given by solid 
red (dashed green) lines. The Fermi energy is at eV. 



E(eV) 
FIG. 6. (color online) Electronic density of states (heavy 
red line) and projections onto the LAPW spheres of Os d 
character (dashed green line) and H s character (dotted blue 
line) for Sr20sH6 with the relaxed LDA H positions, on a per 
formula unit both spins basis. 
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FIG. 8. Electronic density of states (solid) and Fe d projec- 
tion onto the LAPW sphere (broken, radius 2.1 ao) for FeBe2 
on a per per formula unit basis. Spin-up is above the and 
spin down is shown as negative. The Fermi level is at eV. 
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